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Abstract — The objective of this paper is to investigate the dynamic characteristics of a flat plate and a beam subjected to a 
moving load with the effects of inertia force, Coriolis force and centrifugal force considered. To this end, the theory of 
moving mass element for plate and that for beam are presented, where the property matrices of the last elements are derived 
based on the superposition principle and the definition of shape functions. It is found that the order of the property matrices 
of the moving mass element for plate is 24x24, vv/zi/e that for beam is 12x12 Combination of the property matrices of the 
moving mass element and the overall property matrices of the plate (or beam) itself gives the overall property matrices of the 
entire structural system. Because the property matrices of the moving mass element have close relationship with the 
instantaneous position of the moving load, they are time-variant and so are the overall property matrices of the entire 
structural system. For validation, the vibration characteristics of the rectangular plate due to a moving load are compared 
with those of the beam, with its sizes being the same as those of the plate, due to the same loading conditions and satisfactory 
agreement is achieved. Some factors closely relating to the title problem, such as the moving-load speed, acceleration, 
inertia force, Coriolis force and centrifugal force, are investigated. Numerical results reveal that all the above-mentioned 
parameters affecting the dynamic responses of the plate to some degree. 
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I. INTRODUCTION 

Flexible structures undergoing external loadings are important research topics. For the structures with slender shapes, the 
theories of beams are usually adopted in the analyses. However, if the structures have plate -like (rather than beam-like) 
shapes, the theories of beams are no longer available for accurately estimating the dynamic behaviour of the structures. For 
this reason, many researchers have studied the vibration characteristics of structures by means of the plate theories [1-14]. 
For instance, Manoach [1] has studied the dynamic behaviour of elastoplastic thick circular plates due to different types of 
pulses by using the Mindlin plate theory. Wu, Lee and Lai [2] have performed the forced vibration analyses of a flat plat 
under various moving loads with finite element method incorporated with the Newmark direct integration method. Gbadeyan 
and Oni [3], Fryba [4], Lin [5], Shadnam, Rofooei and Mehri [6] and Renard and Taazount [7] have investigated the dynamic 
behaviour of beams and plates subjected to moving forces and moving masses. Marchesiello et al. [8] and Chatterjee and 
Datta [9] have, respectively, studied the dynamic behaviour of multi-span bridges and arch bridges under moving vehicle 
loads, where the bridges are modelled as plates. Takabatake [10] has presented a simplified analytical method for calculating 
the dynamic responses of a rectangular plate with stepped thickness and subjected to moving loads. Rossi, Gutierrez and 
Laura [11] have studied the forced vibration responses of a rectangular plate undergoing a stationary distributed harmonic 
loading. Shadnam, Mofid and Akin [12] have formulated the forced vibration problem of a rectangular plate due to a single 
force (or mass) moving along an arbitrary trajectory by means of the analytical and numerical approaches. Wu [13] and 
Kononov and Borst [14] have researched the vibration characteristics of a plate due to forces moving along a circular path, 
respectively, using the finite element method and the analytical approach. From the literature listed above, it is found that the 
effects of inertia force, Coriolis force and centrifugal force, induced by the moving load, cannot easily take into account the 
dynamic responses of the plate. To solve this problem, the theory of moving mass element is presented in this paper. 

Firstly, under the assumption that the moving load is regarded as a concentrated mass, the property matrices of the moving 
mass element for plate are derived based on the superposition principle and the definition of shape functions. Because the 
property matrices of the moving mass element have something to do with the instantaneous position of the moving load, the 
last matrices vary with time. Adding the property matrices of the moving mass element to the overall property matrices of the 
plate itself yields the time-variant overall property matrices of the entire structural system. For validation, the property 
matrices of the moving mass element for beam are also derived and the vibration characteristics of the rectangular plate due 
to a moving load are compared with those of the beam, with its sizes being the same as those of the plate, due to the same 
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loading conditions and satisfactory agreement is achieved. Some pertinent factors closely relating to the title problem, such 
as the moving-load speed, acceleration, inertia force, Coriolis force and centrifugal force, are studied. Numerical results 
reveal that all the above-mentioned parameters have significant influences on the dynamic responses of the plate. 

II. Property matrices of moving mass element for plate 

Figure 1 shows the s'* plate element subjected to a moving concentrated mass m c at the instant of time t. Since the 
concentrated mass m c is moving along a vibrating path on the plate, the vertical (z) velocity and acceleration of the moving 
mass are respectively given by 
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where w- = w- (x, y, t) , w-(x, y,t) and w-{x,y,t), respectively, represent the vertical ( z ) displacement, velocity and 
acceleration of the moving mass m c at position ( x , v ) and time t ; x and y , respectively, represent the velocities of the 
moving mass m c in the x and y directions (cf. Figure 2); while x and y , respectively, represent the accelerations of that 
in the x and y directions. It is noted that the mass is assumed to be in close contact with the plate at all time. 





Figure 1 The s* plate element subjected to a moving load at the instant of time t . 


For convenience. Equation ( lb) is re-written as 
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where the superscripts x , y and dot, respectively, represent the derivatives with respect to x , y and time t ; while V mx = x , 
v n, y = y . y„ = x and V my = y . 
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If the plate is vibrating, then x , y and z force components of the contact point, between the moving mass and the plate 
element, induced by the vibration and curvature of the plate element are, respectively, given by 


/*(*>: y>t) = m c w-8(x-V mx t)S(y-V my t) (3a) 

f-y (*> y, t ) = m c w- y 8(x - V mx t )8( y - V my t) (3b) 


f-M y,t) = m c [wfVl + 2wfV mx V my + 2wy m 


+ 2w2V my+W fV^+w^+wlV my +w- z ]S(x-V mx t)S(y-V my t) 


(3c) 


where S(x-V mx t) and 8(y-V my t) represent the Dirac delta function; while w- = w-(x,y,t) and w y = w y (x, y,t ) , 
respectively, represent the x and y accelerations of the plate element at position ( x , y ) and time 1 . 

In such a case, the equivalent nodal forces, f x p) ( k =1 to 24), of the s' p plate element are given by [13] 
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where <t> ( k p> ( k =1 to 24) are shape functions with the non-zero ones given by [13] 
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where i px and t are respectively the length and width of the rectangular plate element (see Figure 1), whereas x m p (t) and 
y ( „ p, (t) are respectively the local x and y positions of the concentrated mass m c at time t. 


Based on the superposition principle and the definition of shape functions, the local x , y and z displacements of contact 
point, vv x , w and w , , can be obtained from 
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where u (p) ( i =1 to 5, 7 to 1 1, 13 to 17 and 19 to 23) are the nodal displacements of the nodes of the plate element on which 
the moving concentrated mass m c applies. 

Introducing Equations (3) and (6) into Equation (4), and writing the resulting expressions in matrix form yields 
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In Equation (8), all the coefficients of [m (1) ] 24x24 , [c (1) ] 24x24 , [c (2) ] 24x24 , [k (1) ] 24x24 , [& ,2) ] 24x24 , [& <3) ] 24x24 , [k {4) ] 24x24 and 
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It is noted that the superscripts x and y represent the derivatives with respect to x and y , respectively. 

In Equation (8), [m (p) ] , [c tp) ] and [ k {p) ] are the mass, damping and stiffness matrices of the presented moving mass element 
for plate. Since the last matrices have something to do with the shape functions of the plate element on which the moving 
mass applies (see Equations (5) and (9)), their properties change from time to time if the position of concentrated mass m c 
on the plate is time-variant. 

III. Dynamic responses of a flat plate due to a moving load 

For a multiple degree of freedom damped structural system, its equations of motion is given by [15] 

[M ] { W) } + [C ] {£(f ) } + [K] {q{t) } = {F (f) } (10) 

where [M], [C] and [A] are, respectively, the overall mass, damping and stiffness matrices; {g(f)}, {q(t)} and {<?(0} 
are, respectively, the acceleration, velocity and displacement vectors, whereas [F (t)) is the overall external force vector at 
any time t. For a plate subjected to a moving load, the overall property matrices, [M] , [C] and [A] , and the external force 
vector {F(t)} are determined using the expressions introduced in the subsequent subsections. 


3.1 Overall property matrices of the entire structural system 

To take the effects of inertial force and centrifugal force of the moving load into account, the overall mass matrix, [ M ] , and 
stiffness matrix, [A], of the entire structural system must be determined by adding the mass and stiffness matrices, [m p) ] 
and [fc (p) ] , of the moving mass element to the overall ones of the plate itself, i.e., 

[M] = [M p ] + [m (p) ] 24x24 (1 1a) 

[K] = [K p ] + [k (p) ] 2 4x24 (lib) 

and imposing the prescribed boundary conditions. In which, [M ] and [K p ~\ are, respectively, the overall mass and stiffness 
matrices of the plate itself and may be obtained by assembling all its element mass and stiffness matrices [16]. 

In Equation (11), all the coefficients of the mass and stiffness matrices, [M] and [A], are exactly the same as the 
corresponding ones of [M ] and [^ K ] , i.e., 

M,j=M pJj (12a) 

K ij= K P .ij (1,7 = 1 ton) (12b) 


except 


M D:Dj =M pADj +m< p) 

K DlDj =K p DiDj +k^ (1,7=1 to 24) 


(13a) 

(13b) 


In Equations ( 12)-( 1 3), n and IX (1=1 to 24 ) are, respectively, the total degrees of freedom of the entire structural system 

and the numberings for the 24 degrees of freedom for the four nodes of the D' h plate element on which the moving load 
applies at time t. 

Due to the fact that the elementary damping matrix of the structural system is difficult to find from the existing literature, the 
overall damping matrix [C ] of the plate itself is assumed to be proportional and determined by using the Rayleigh damping 
theory [15]. 
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where £ and represent the damping ratios corresponding to any two natural frequencies of the structural system, 0 ) and 
co j , respectively. 

In such a case, the overall damping matrix [C] of the entire structural system can be obtained from 
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except 
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where the meanings of n and I) l ( 1 = 1 to 24 ) are exactly the same as those of Equations (12) and (13). 


Because the mass, damping and stiffness matrices, [m <p> ] , [c <p) ] and [£ (p) ] , of the moving mass element for plate are time- 
dependent, so are the overall mass, damping and stiffness matrices, [M ] , [C] and [K ] , of the entire structural system. 

3.2 Overall external force vector 

Figure 2 shows a concentrated mass m c moves, with acceleration a , on a rectangular plate. If, at any instant of time t, the 
concentrated mass locates at the position 3c ,„(t) and y„(t) of the plate, then the external forces on all the nodes of the 
plate are equal to zero except the four nodes of the s'* plate element at which the concentrated mass is located at time t. 


where 
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In Equation (16), f s p (t) (1=1 to 24) are the equivalent nodal forces on the 24 degrees of freedom of the j ' plate element. 
In Equation (17), Equations (17a) and (17b) are those due to the horizontal (J) and horizontal (y) inertia forces of the 
moving load, respectively; while Equation (17c) are those due to the external load m c g induced by the concentrated mass 
m c and located at the local coordinates [ x , y ]=[ x <p) (t ) , y { , p) (t) ] of the s'* plate element. It is noted that the symbols g 
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and 9 appearing in the last expressions are respectively the gravity acceleration and the angle between the direction of 
moving concentrated mass m c and the x axis (see Figure 2). 



Figure 2 A concentrated mass " lc moves, with acceleration a , on a flat plate. 

3.3 Solution of equations of motion 

In the present study, the dynamic responses of the flat plate due to a moving load, with the effects of inertial force, Coriolis 
force and centrifugal force considered, are determined by solving the equations of motion of the entire structural system 
given by Equation (10). The solution procedures are described in the following. 

1. Since the initial conditions of the structural system are assumed to be “at rest” in this paper, [g(0)} = {0} , {^(0)} = {0} 
and (|(0)} = {0}. 

2. Using Equations (8a)-(8c) to calculate the property matrices of the moving mass element for plate, [m (p) ], [c ,p) ] and 

[k ' p) ~\ , at any time t i+1 = t i + At ( i =0,1,2,...). Where Ar is time interval used for calculating the dynamic responses of 

the structure and taken to be 0.001s in this paper. 

3. Using Equations (11)-(13), and imposing the prescribed boundary conditions, to determine the overall mass and 
stiffness matrices of the entire structural system, [M] and [A] , at any time t i+l = r ; + Ar (i =0,1,2,...). 

4. Determine the natural frequencies, co t ( i =1,2,. . .), of the entire structural system by means of Jacobi algorithm [16]. 

5. Perform the Rayleigh damping theory (see Equations (14a)-(14c)) to determine the overall damping matrix, [C ], of 

the plate itself and calculate the overall damping matrix, [C ] , of the entire structural system (see Equation (15)). 

6. Evaluate the overall equivalent nodal force vector (see Equations (16)-(17)) of the structure due to moving load at any 
time t M = tj + At ( i =0,1,2,...). 

7. Solve the equations of motion of the entire structural system. Equation (10), with Newmark direct integration method 
[16] for its dynamic responses at any time t i+l = f, + Ar ( i =0,1,2,...). 

8. Repeat steps 2-7 to obtain the dynamic responses of the structural system at any time t M = r ; + Ar ( i =0,1,2,...). 
Because the property matrices, [rrr <p> ] , [c <p> ] and [k <,,) ] , of the moving mass element for plate are time-variant, so are the 
overall mass, damping and stiffness matrices of the entire structural system. In such a case, the overall mass, damping and 
stiffness matrices must be calculated at each time step (see steps 3-5). For this reason, the time required for the computer 
calculation by the present technique is greater than that required by the existing literature without considering the inertial 
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force, Coriolis force and centrifugal force of the moving load(s). However, it is believed that this should be the cost that one 
should pay if more satisfactory results are hoped to obtain. 

IV. Dynamic responses of the beam due to a moving load 

The dynamic responses of the flat plate due to a moving load can be determined by means of the formulations presented in 
sections 2-3 of this paper. Because the vibration responses of the flat plate subjected to a moving load should be close to 
those of the beam, with its sizes being the same as those of the flat plate, subjected to the same loading conditions [2], the 
formulations for calculating the dynamic responses of the beam due to a moving load are also presented in this section. Thus, 
the vibration characteristics of the rectangular plate due to a moving load obtained from the formulations of sections 2-3 can 
be validated by using those of the beam obtained from the formulations of this section. 



Figure 3 The Sb beam element subjected to a moving load at the instant of time 1 . 


Since each beam element consists of two nodes and 12 degrees of freedom (see Figure 3), the mass, damping and stiffness 
matrices, [m (W ], [c <i) ] and [£ (6) ] , of the moving mass element for beam are 12x12 matrices and can be derived by using 
the similar procedures of section 2, where all the coefficients of the above-mentioned matrices are equal to zero except that 

mf =m c <S>f ) <S>f (i. 

7=1,7) 

(18a) 

m™=r. (i. 

j =2,6,8,12) 

(18b) 

mf = (i. 

7=3,5,9,11) 

(18c) 

cT= 2m c V mx ^T 

( i, j =2,6,8, 12) 

(18d) 

= 2my mx &r&’; ) 

(*, 7 =3,5,9,11) 

(18e) 


(i,;=2,6,8,12) 

(18f) 


Page | 8 




International Journal of Engineering Research & Science (IJOER) 


[Vol- 1 , Issue-8, November- 2015] 


C = ( i, J =3,5,9,11) (18g) 

where V mx is the velocity of the moving concentrated mass m c in the 3c direction; while <t>f (1) (1 =1 to 4) are the shape 
functions with non-zero ones given by [15] 

<*>f = !-£», &?=&? = l-3g„ 2 + 2g b \ Of =$»}« = 0, 

of (t) = of (0 = [ ?6 - 2 ? f + ]g „ , Of = g b , 

Of = Of = 3f 2 - 2gf , <f>f = <t>f = + g b 3 Vb (19a) 

s t =x?(f)/e b d9b) 

where xf (f ) is the distance between the position of the moving mass m c and the left-end of the beam element on which it 
applies; while l b is the length of the beam element (see Figures 3 and 4). 

To determine the overall property matrices of the entire structural system, one requires to add the property matrices of the 
moving mass element for beam and those of the entire beam structure together, i.e., 


[M] = [MJ + [m (i) ] 1M2 (20a) 

[C]=[CJ + [c w ] 12xl2 (20b) 

[K] = [K h ] + [k ,b) l llx]2 (20c) 

where 

M ij = M b.a ( h j = 1 10 n ) (2 1 a ) 

C„=C hll (f, ; = 1 ton) (21b) 

K, K hii (1,7 = 1 ton) (21c) 

except that 

M DiDj = M bADj + mf (*,7=1 to 12) (22a) 

C m =C b M>,+ c f ( i,7=l to 12) (22b) 

k d,d, = K b,D,D, + m T ( i,j = 1 to 12 ) (22c) 


In the last expressions, [M b ] , [C b ] and [If] are the overall mass, damping and stiffness matrices of the beam itself, 
respectively. In which, the overall damping matrix, [ C b ], is obtained by means of the Rayleigh damping theory [15] (cf. 
Equation (14)). Moreover, D (1=1 to 12) respectively represent the numberings for the 12 degrees of freedom of the two 
nodes of the D' h beam element on which the moving load applies at time t. 

If, at any instant of time t, the concentrated mass m c is located at the position xf (t) of the beam (see Figure 4), then the 
overall external force vector of the entire beam induced by the concentrated mass m c takes the form 

{F(t)} = [0 0 - f?\t) f£\t) •••/;>) O) ■■■ 0 Of (23) 

where 

/f (t) = m c g^f ) (1=3, 5, 9, 11) (24a) 
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f ( s b> (t) = 0 (i =1, 2, 4, 6, 7, 8, 10, 12) (24b) 

It is worthy of mentioned that the symbol g refers to the acceleration of gravity and the subscripts, s, (i= 1 to 12) in 
Equations (23)-(24) represent the 12 degrees of freedom of the s b beam element at which the concentrated mass m c is 
located. 



Figure 4 A beam subjected to a moving load. 

Finally, one may use the similar procedures present in section 3.3 to solve the equations of motion. Equation (10), of the 
structural system for the vibration responses of the beam due to a moving load. 

V. Numerical results and discussions 

In this section, one plate model and one beam model, respectively, subjected to a moving load are studied. In which, the plate 
is made of steel with mass density p p = 7820 kg/m 3 , modulus of elasticity E p = 206.8 GN/m 2 and Poisson’s ratio V = 0.29, 

while its dimensions are: length L px =1.0 m, width L = 0.5 m and thickness L p = 0.01 m (see Figure 5). The entire plate 

is modelled with 32 identical 0.125 mx 0.125 m rectangular plate elements and 45 nodes. For conveniences, the last 
rectangular plate with side AB and side CD being constrained as pin joints is called Pin-plate. 



Figure 5 A pinned-pinned rectangular plate (Pin-plate) subjected to a moving load with mass 

= 2.3 KG AND A CONSTANT SPEED V =10.0 M/S. 

Since the vibration characteristics of the rectangular plate, with its bending effects neglected (i.e., Poisson’s ratio V =0), will 
be close to those of the beam [2], a beam model is also studied in order to validate the availability of the presented theory. 
The beam model corresponding to Pin-plate is called Pin-beam, where the dimensions and material properties of Pin-beam 
are exactly the same as those of Pin-plate (i.e., length L b = 1.0 m and cross-sectional area 0.5 mx 0.01 m, mass density p h = 

7820 kg/m 3 , modulus of elasticity E b = 206.8 GN/m 2 ). In addition, the finite element model for Pin-beam is composed of 9 
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nodes and 8 beam elements (see Figure 6). Three values of Poisson’s ratio ( V = 0.0, 0.15, 0.29) for Pin-plate were studied, 
however, according to the beam theory, the Poisson’s ratio for Pin-beam is equal to zero. Unless specially stated, the 
damping ratios ( ) used for the plate model and the beam model of this paper are taken to be 0.005 and the time interval 
used for calculating the dynamic responses of the entire structural system is 0.001s, i.e.. At = 0.001 s. 



Figure 6 A pinned-pinned beam subjected to a moving load with mass 2.3 kg and a constant 

SPEED V =10.0 M/S. 


5.1 Validation 


The lowest few natural frequencies of Pin-plate (with Poisson’s ratio V =0.0, 0.15 and 0.29, respectively) and Pin-beam are 
listed in columns 2, 3, 4 and 5 of Table 1. Because the total degree of freedom of Pin-plate is much more than that of Pin- 
beam, one can find only four of the lowest ten natural frequencies of Pin-plate to be corresponding to the lowest four natural 
frequencies of the Pin-beam, as shown in the fifth column of Table 1. It is noted that the correspondence between the natural 
frequencies of Pin-plate and those of Pin-beam must be made based on their corresponding mode shapes. From the table, it 
can be found that the natural frequencies of Pin-plate will be close to the corresponding ones of Pin-beam if the Poisson’s 
ratio ( V ) approaches zero. Based on this result, one may infer that the central displacements of Pin-plate subjected to a 
moving load will be close to those of Pin-beam subjected to the same loading conditions if the Poisson’s ratio (v) of the 
plate approaches zero. 


Table 1 

The lowest few natural frequencies oj (Hz) of Pin-plate and Pin-beam. 


Mode No. for Pin- 
Plate 

Pin-plate 

Pin-Beam 

V = 0.29 

V = 0.15 

© 

© 

II 

1 st 

23.6032 

23.4487 

23.3869 

23.3150 

2 nd 

63.1528 

66.1188 

70.0181 


3 rd 

96.3655 

94.9433 

94.4052 

93.2379 

4 th 

146.4804 

150.3048 

155.8149 


5 th 

215.1707 

212.2335 

213.6403 


6 th 

221.5174 

217.402 

215.8878 

209.8451 

rj th 

264.8742 

267.1878 

272.0683 


8 th 

283.4873 

285.2118 

292.2830 


9 th 

380.7570 

385.1473 



10 th 

403.5357 

395.9960 

393.2355 

373.7035 
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If a load with mass m c = 2.3 kg moves with a constant speed V c = 10.0 m/s from side AB to side CD, along the centrelines, 
of Pin-plate (see Figure 5) and Pin-beam (see Figure 6), then the vertical ( z ) central displacements of Pin-plate and Pin- 
beam are shown in Figure 7, where the solid line ( — ■), the solid line with triangles ( — A — ) and the solid line with circles 
( — o — ) represent the vertical ( z ) central displacements of Pin-plate with Poisson’s ratio V = 0.0, 0.15 and 0.29, 
respectively, while the dashed line ( — ) represents those of Pin-beam. From the figure, it is seen that the curves for Pin-plate 
are close to the curve for Pin-beam when the value of Poisson’s ratio (v) approaches zero. This agrees with the numerical 
results presented in reference [2]. 



Figure 7 Vertical ( z ) central displacements of Pin-plate and Pin-beam. 

From all the numerical results presented in this subsection, it is believed that the presented theory is available for the title 
problem and will be used for further studies in this research. 

5.2 Influence of moving-load speed 



Figure 8 A pinned-pinned plate (Pin-plate) subjected to a moving load with mass '” c 2.3 kg and 

v 


MOVING SPEED 
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In this subsection, a load of mass m c =2.3 kg moves with a constant speed V from point E ( J, =0.2 m, y, =0. 1 2m) to point F 
( 3c 2 =0.8 m, y 2 =0Am) of the Pin-plate (see Figure 8) is investigated. Figure 9 shows the time histories for the vertical (z) 
displacements of the centre point of Pin-plate, where the solid curves with circles ( — O — ■) represent the time histories with 
moving-load speed V = 5.0 m/s, those with crosses ( — + — ■) represent the ones with V = 10.0 m/s and those with triangles 
( — A — ) represent the ones with V = 20.0 nils. From the figure, one sees that the larger the moving-load speed, the larger the 
maximum vertical ( z ) central displacements of the flat plate. 

2.00E-005 
1 .00E-005 
0.00E+000 

£ 

o -1 .00E-005 

Q. 

| -2.00E-005 

0 

1 -3.00E-005 

CD 

E 

S -4.00E-005 

Q. 

(/> 

75 -5.00E-005 

c 
CD 

^ -6.00E-005 

"cC 

1 -7.00E-005 

CD 
> 

-8.00E-005 
-9.00E-005 
-1 .00E-004 

0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00 

VV,- y -y f / Vf >,->,)■’+( y,-y,r 

Figure 9 Vertical ( z ) central displacements of Pin-plate due to a moving load with a 
constant speed: (a) v = 5.0 m/s, (b) v =10.0 m/s and (c) V =20.0 M/S. 

5.3 Influence of acceleration 

If the mass m c moves, with initial velocity V =0.0 mis and constant accelerations a =1.0, 4.0 and 8.0 m/s 2 , from point E 
(3Cj=0.2 m, y 1 =0A2m ) to point F (3c 2 =0.8 m, y 2 =QAm) of the Pin-plate, then the time histories for the vertical (z) 
displacements of the centre point of Pin-plate are shown in Figures 10(a), 10(b) and 10(c), respectively. In which, the solid 
curves with circles ( — O — ■), the solid curve with crosses ( — + — ) and the solid curves with triangles ( — A — ) represent the 
responses corresponding to the accelerations a =1.0, 4.0 and 8.0 m/s 2 , respectively. From the figures, one finds that the larger 
the acceleration of the moving mass m c , the larger the maximum vertical (z) displacement of the centre point of Pin-plate. 
Since the acceleration of the moving mass has a close relation with its velocity and the latter significantly affects the vertical 
( z ) responses of the Pin-plate, one must also consider the effect of velocity in addition to the acceleration. 

5.4 Influence of inertia force 

From the formulation of this paper, one can see that the effect of inertia force induced by the moving mass is to appear in the 
mass matrix [m lp) ] of the moving mass element for plate as one may see from Equation (8a). Hence, if the mass matrix of 
the moving mass element is taken to be zero, i.e., [m ip> ] = [0], then the effect of inertia force due to the moving mass will 
disappear. 
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Figure 10 Influence of acceleration on the vertical ( z ) central displacements of Pin-plate 
DUE TO A MOVING LOAD WITH A CONSTANT ACCELERATION: (A) a = 2.0 m> ^ , (B) a =4.0 m/s 2 AND (C) 

a = 8.0 mls \ 
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(A) 



(B) 



(C) 


Figure 11 Influence of inertia force on the vertical ( z ) central displacements of Pin-plate 
DUE TO A MOVING LOAD WITH A CONSTANT SPEED: (A) V = 5.0 M/S, (B) V =10.0 M/S AND (C) V =20.0 M/S. 
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Figure 12 Influence of Coriolis force on the vertical ( z ) central displacements of Pin-plate 
DUE TO A MOVING LOAD WITH A CONSTANT SPEED: (A) V = 5.0 M/S, (B) V = 10.0 M/S AND (C) V = 20.0 M/S. 





International Journal of Engineering Research & Science (IJOER) 


[Vol-1, Issue- 8, November- 2015] 


2.00E-005 


3s o.ooe+ooo m 


$ -2.00E-005 


co -4.00E-005 


-6.00E-005 


« -8.00E-005 


-1.00E-004 


with centrifugal force considered 

-£) — with centrifugal force neglected 




o 




/ \ 


6 


k 


© 


o 


o 

r 


\ / 
V 


\ ,/ 

kj 


0.00 0.02 0.04 0.06 0.08 0.10 0.12 

Time (sec) 


(A) 



(B) 



(C) 

Figure 13 Influence of centrifugal force on the vertical ( z ) central displacements of Pin- 
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Figure 1 1 shows the time histories for the vertical ( z ) displacements of the centre point of Pin-plate when the mass m c =2.3 
kg moves with a constant speed V from point E ( 5j =0.2 m, y, =0. 12m) to point F ( x 2 =0.8 m, y 2 =0.4 m) of the flat plate. In 
which, the curves with circles ( — O — and — O— ), crosses ( — + — and and triangles ( — A — and —A—) represent the 

time histories with moving-load speed V =5.0 nils, V = 10.0 mis and V =20.0 mis, respectively. Besides, the solid and dashed 
curves represent those with the effect of inertia force of moving load considered and neglected. From the figure, one sees that 
the maximum vertical ( z ) central displacement with the effect of inertia force considered is larger than that with the effect of 
inertia force neglected. Thus, the effect of inertia force is important and should be considered in the formulations. 

5.5 Influence of Coriolis force 

Similarly, one can ignore the effect of the Coriolis force due to moving load by taking the damping matrix of the moving 
mass element for plate to be zero, i.e., [c (p> ] = [0], The same example as that of the last subsection is studied and the vertical 
(z) displacements of the centre point of the flat plate are shown in Figure 12. The legends for the curves in the figure are 
exactly the same as those in Figure 1 1 except that the inertia force is replaced by the Coriolis force. From the figure, one sees 
that the Coriolis force affects the vertical ( z ) central displacements of the plate to some degree. 

5.6 Influence of centrifugal force 

In this subsection, the effect of centrifugal force due to moving load is ignored by taking the stiffness matrix of the moving 
mass element for plate to be zero, i.e., [k <p) ] = 1 0] . The same plate as that of the last subsection is investigated and the 
vertical ( z ) and central displacements of the Pin-plate are shown in Figure 13. The legends for the curves are exactly the 
same as those in Figure 12 except that the Coriolis force is replaced by the centrifugal force. 

From Figure 13, one sees that the influence of the centrifugal force on the vertical (z) central displacements of the plate 
increases with increasing the moving-load speed. This is because the magnitude of the centrifugal force appearing in the 
stiffness matrix [k (,,) ] of the moving mass element for plate is proportional to the square of the moving-load speed (see 
Equation (8c)). 


VI. Conclusion 

1. To take account of the effects of inertia force, Coriolis force and centrifugal force of the moving load, the theory of the 
moving mass element for plate and that for beam are presented. In which, the property matrices of the moving mass 
elements are derived based on the superposition principle and the definition of shape functions. It is found that the order 
of the property matrices of the moving mass element for plate is 24x24 , while that for beam is 12x12 . 

2. Combination of the property matrices of the moving mass element and the overall property matrices of the plate (or 
beam) itself gives the overall property matrices of the entire structural system. Because the property matrices of the 
moving mass element for plate and that for beam have something to do with the instantaneous position of the moving 
load, both the property matrices of the moving mass element and the entire structural system are time-variant. 

3. The moving speed, acceleration, inertia force, Coriolis force and centrifugal force of the moving load have significant 
influences on the vertical ( z ) dynamic responses of the flat plate. Thus, all the above-mentioned parameters should be 
considered in the formulations. 
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